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ABSTRACT 

The  near  field  pressure  and  the  radiation  impedance  for 
an  infinite  phased  array  of  circular  pistons  (with  the  same  velo- 
city magnitude)  in  a rigid  plane  baffle  are  calculated.  A slowly 
convergent  infinite  series  expression  for  the  pressure,  obtained 
from  the  appropriate  Green's  function,  is  transformed  by  the  use 
of  the  Poisson's  sum  formula  into  a more  rapidly  convergent  in- 
finite series  expression.  Numerical  results  are  presented.  A 
rough  estimate  of  the  location  and  value  of  the  maximum  pressure 
in  the  near  field  is  obtained  under  certain  restrictions.  The 
radiation  impedance  of  a piston  in  the  infinite  array  agrees 
well  with  the  average  radiation  impedance  of  a piston  in  a large 
finite  array. 
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A.  INTRODUCTION 


A sonar  designer  is  frequently  faced  with  the  task  of 
calculating  or  estimating  the  peak  pressure  in  the  near  field  of 
a phased  array  because  the  usable  power  output  of  an  acoustic 
underwater  transducer  is  limited  by  the  requirement  that  the 
acoustic  pressure  in  water  should  not  exceed  the  value  at  which 
cavitation  begins.  The  computation  of  the  peak  pressure  can  be 
a rather  formidable  job  because  the  summation  of  the  pressures 
due  to  each  array  element  usually  cannot  be  reduced  to  anything 
simple  analytically,  and  the  location  of  the  maximum  pressure  is 
not  known  a priori.  Furthermore,  to  choose  the  optimum  sonar 
design  one  would  like  to  know  how  the  pressure  maximum  changes 
when  the  sonar  array  parameters  (size  of  elements,  spacing  bet- 
ween elements,  number  of  elements,  direction  to  which  the  array 
is  steered  in  the  far  field,  etc.)  are  varied,  which  might  mean 
that  the  computation  of  the  near  field  pressure  has  to  be  per- 
formed for  a large  number  of  arrays.  Thus  any  simplifications 
are  welcome  as  long  as  the  simplifications  permit  to  study  the 
effects  of  some  of  the  sonar  parameters  on  the  near  field 
pressure . 

Several  calculations  of  the  near  field  pressures  for 
special  distributions  of  elements  in  arrays  are  available.* 

For  an  infinite  array  two  advantages  in  the  representation  of  the 
near  field  pressure  are  attained:*^  1)  one  has  to  consider  the 
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pressure  in  front  of  one  element  only,  because  due  to  symmetry 
all  elements  are  the  same,  except  for  a phase  factor;  2)  the 
infinite  sum  of  pressures  due  to  each  array  element  can  be  trans- 
formed by  the  use  of  the  Poisson's  sum  formula  into  a more  ra- 
pidly convergent  series  of  terms  which  resemble  waveguide  mode 
terms  in  the  theory  of  electromagnetic  waves. 

One  hopes  that  the  inferences  drawn  from  the  infinite 
array  calculations  about  the  effects  on  the  near  field  due  to 
changes  in  the  element  sizes  or  spacings  between  elements  will 
also  apply  to  arrays  with  a finite  number  of  elements.  The  two- 
dimensional  problem  of  an  infinite  array  of  infinite  vibrating 
strips  in  a rigid  plane  baffle  has  been  studied  before  we 
will  now  consider  the  three-dimensional  problem  of  an  infinite 
phased  array  of  circular  pistons  in  a rigid  plane  baffle.  The 
radiation  impedances  of  the  pistons  will  also  be  calculated.  The 
velocity  magnitudes  of  the  pistons  are  assumed  all  to  be  equal, 
but  the  phases  are  adjusted  so  as  to  steer  the  far  field  pressure 
maximum  in  a specified  direction. 

B.  THE  PRESSURE  FIELD 

Let  the  infinite  array  and  baffle  occupy  the  plane 
z = 0,  (see  Fig.  1).  The  piston  centers  are  located  at  x = nd^ 
y = md 2,  n = 0,*1, ±2,...;  m =»  • The  velocity  (in 

the  z-direction)  of  the  piston  in  the  mtl1  row  and  nC^  column  is 
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assumed  to  be 


v 

mn 


loot  + iy 


ve 


mn 


(1) 


where  v is  a complex  constant,  cjd  the  radian  frequency,  t the  time, 
and 


7,,^  = -mkd2sin0osin0o  - nkd1sin0Qcos0o , 


(2) 


where  k = co/c,  c is  the  velocity  of  sound  in  water,  and  0Q, 

0Q  are  angles  in  spherical  coordinates  which  determine  the  di- 
rection to  which  the  far  field  pressure  maximum  is  steered. 

If  the  actual  pressure  at  the  observation  point  due  to 
th  th 

the  m , n piston  alone  is  p , let  P be  a dimensionless 

mn  mn 

pressure , 


Pmn  - (3) 

where  p is  the  density  of  water.  Similarly  let  P be  the  tota.' 
dimensionless  pressure  at  the  observation  point  due  to  all  pis- 
tons , 


00 

-I 

nfc=-oo 


Pe 
mn 


*7, 


mn 


3 


(4) 
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By  the  use  of  a Green’s  function  can  be  written 


.11,12 


2tt  a 


*-  n « _ ilcR.  I 

Pmn  = (ik/27r)  X da°  Jodr°  ^ 6 'R 


th  th 


where  R is  the  distance  from  a point  on  the  m , n piston 
mn 

at  r , a to  the  observation  point  at  x,  y,  z;  see  Fig.  lj 
o o 


>W  ■ <»2  + >1  ♦ 


xn  = X - ndi  - rocosao) 


y = y - md0  - r sina . 
•ym  J 2 o o 


We  now  use  the  following  integral  representation  for  spherical 


waves 


-ikR_ 


OO  00 


\nn  2tt 


I ds-,  I ds2  / ds. 

J.oo  ^J—00  ^-00 


-islxn‘is2ym'is3z 


S-J^  4"  S 2 4"  ■ k 


which,  when  integrated  over  s3>  becomes 
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- ikR 


mn 


X 


mn 


■ kf  *1 f 

^—oo  ■'-a 


ds, 


-tslxn-is2ym-qz 


(10) 


where 


q = < 


i(k2  - s2 
(s2  + s2 


s2)1/2*  0 < si  + s2  1 k2> 
k2)1/2,  s2  + s2  > k2. 


(11) 


If  we  combine  Eqs.  (4),  (5),  and  (10),  and  interchange 
summation  and  integration,  we  obtain 


2tt  a 


. , Z.TT  a 00  00  ' oo  oo 

P / d°°/dr°ro  Z IJdSl/dS2 

-/Q  ^O  m=  -00  n=_oo-'-oo  -oo 


(l/q)exp(i7mn  - iS;Lxn  - is2ym  - qz) 


(12) 


We  now  apply  Poisson's  sum  formula*^’ ^ 


OP  00 

£ f ds  f(s)e 

n=  - oo  oo 


2frins 


OO 

I f<"> 


n=-oo 


(13) 


twice  to  Eqs.  (12),  which  yields 
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I 

I 

I 


I 


I 

I 

I 


i 


. “ °°  -T(n,m)kz  2tt  a 

P - E t ' ' T(n»  f d“o/droro 

n=-°°  nt=-°°  Jo  Jo 

exp ^ih(n) (krocosaQ  - kx)  + ig(m) (krosinaQ  - ky)]  . 


h(n)  = sin0Qcos0  + 2Tm/kd^ 


(14) 

(15) 


g(m)  = sin0osin0o  + 27rm/kd2 


w(n,m)  = |[h(n)]  2 + [g(m)]  2^2 


T(n,m)  = 


ijl-  [wCn.m)]2}1/2  , 0 < w(n,m)  < 1; 
1/  [w(n,m)J2  - l}1^2,  w(n,m)  >1. 


Since 


16 


J da0  exp|ikro[h(n) 


cosaD  + g(m)sinaOJ 


= 2tt  Jo[w(n,m)kro]  , 


(16) 

(17) 


(18) 


(19) 
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where  J is  a Bessel  function  of  order  zero,  and  since 
o ’ 


17 


a 

f droro  Vwkro>  = (a/wk)Ji(wka) 

Jo 


(20) 


where  is  a Bessel  function  of  order  one,  we  finally  obtain 
from  Fq . (14) 


2 00  00 


p = 


i-jra  r*  ^llw^ 

~ L L»  w(n, 


2Ji  £w(n,m)kaj  e-ih(n)kx-ig(m)ky-T(n,m)kz 


m)ka 


T(n,m) 


n=-°°  m=-°° 


(21) 


which  is  the  desired  exact  expression  for  the  pressure  at  x,y,z 
(either  near  field  or  far  field)  due  to  the  infinite  array. 

The  form  of  Eq . (21)  is  similar  to  the  expansions  of 

18 

electromagnetic  fields  in  waveguide  modes.  The  double  series 
in  Eq . (21)  is  more  rapidly  convergent  than  the  original  series 
from  Eqs.  (4)  and  (5);  for  z > 0 the  terms  in  the  sum  are  de- 
creasing exponentially  with  increasing  m and  n,  and  for  z = 0 
the  terms  are  decreasing  as  [(n/dp2  + (m/d2)  2]5//4. 

Let  us  now  consider  some  special  cases  of  Eq . (21). 

If  kz  » 1 and  kd^  < tt,  kd2  < tt,  all  the  terms  are  exponentially 
decreasing  except  the  one  for  n = m = 0.  Since 


w(0,0)  = sin0Q 


(22) 
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and 


T(0,0)  = icos0o> 


(23) 


if  we  neglect  the  exponentially  attenuated  terms,  we  obtain  a 
single-term  approximation  for  the  pressure 


2J1(kasln9Q) 
1 * kasinQ 


Pc,  kz  » 1, 


(24) 


where  Pc  is  the  dimensionless  pressure  which  one  would  have  in 
the  water  if  the  whole  z = 0 plane  were  vibrating  continuously 
with  the  normal  velocity 


io)t-ikx  sin0^  cos0_-iky  sin0  sin0  : 
vc  = ve  ° ° ' ° ° (25) 

i.e.,  if  instead  of  a plane  partially  filled  with  pistons  (each  of 
which  move  with  the  same  velocity  as  a rigid  body)  we  would  have 
a plane  on  which  each  infinitesimal  area  element  moves  with  the 
velocity  given  by  Eq . (25),  then  the  pressure  would  be*^ 

pc  = (l/cos0o)exp(-ikx  sin0ocos0o  - iky  sin0osin0o  - ikz  cos9o) 

(26) 
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Thus  the  approximate  as  compared  to  Pc  for  kz  » 1 is  re- 

2 

duced  by  the  ratio  of  the  piston  area  ira  to  the  area  of  the 

2 

array  per  piston  dj^  (because  only  the  fraction  ~a  /dj^  ° f the 
plane  z = 0 is  vibrating  in  the  piston  array,  while  the  whole 
plane  is  vibrating  for  the  continuous  plane  radiator)  and  by  the 
directivity  of  the  circular  piston  2j^(kasin0Q) / (kasin0Q)  < 1 
(because  the  piston  vibrates  as  a rigid  body,  therefore  in  a 
direction  0Q  4 0 the  wavelets  arriving  from  different  portions 
of  the  same  piston  are  not  quite  in  phase,  which  results  in  par- 
tial cancellation;  all  the  waves  are  in  phase  for  the  continuous 
plane  radiator) . 

The  n = m = 0 term  dominates, and  we  obtain  P^  from 
Eq . (24)  as  an  approximation  for  P also  when  d^  ->  0,  d2  -»  0, 
a ->  0,  i.e.,  when  the  piston  array  consists  of  small  pistons 
close  together,  because  such  an  array  approximates  the  continuous 
plane  radiator. 

Fig.  2 shows  |p|  vs.  kz  for  equal  spacings  between  rows 
and  columns,  kd^  = kd2  = 3,  and  three  values  of  the  piston  radii, 
ka  = 0.5,  1.0,  and  1.5;  the  pressure  is  evaluated  on  the  axis 
through  the  center  of  a piston,  kx  = ky  = 0,  and  the  array  radiates 
in  the  broadside  direction  normal  to  the  plane  containing  the 
array,  0Q  = = 0°.  For  ka  = 1.5  the  pistons  are  touching  each 

other.  For  large  values  of  kz  |p|  approaches  the  constant  |PjJ 
given  by  Eq . (24). 
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Fig.  3 shows  |p  | vs.  kx  for  ©o  = 0°,  60°,  and  80°, 
and  0O  = 0°;  the  other  parameters  are:  kd^  = kd2  = 3,  ka  = 1, 
ky  = kz  = 0.  For  0Q  = 0°  the  pressure  distribution  is  relatively 
smooth,  while  for  larger  angles  it  becomes  more  undulatory.  As 
0Q  increases,  the  maximum  |p|  increases  because  the  leading  n = 
m = 0 term  in  Eq . (21)  is  proportional  to  1/T(0,0)  = l/icos0Q. 

For  0Q  ->  90°  we  would  have  |p|  -»  °°. 


The  single-term  approximation  P^  does  not  show  any 
undulatory  behavior;  therefore  let  us  obtain  a two-term  approxi- 
mation P2  to  explain  the  maxima  and  minima  in  Fig.  3.  For 
simplicity  let  0Q  = 0°,  kd^  < 7 r,  kd2  < kd^,  and  0q  >0°,  then 
the  next  largest  term  after  the  n = m = 0 term  is  the  n=-l, 
m=0  term.  Let  us  call  the  approximation  consisting  of  the  two 
terms  n=-l,  0 and  nt=0  the  pressure  P2.  Let  us  further  assume 
that  ka  is  small  enough  so  that  2J^  £w(n ,m) kaj/w(n ,m) ka  «1  for 
both  terms.  Thus  we  obtain  an  approximation  which  for  0O  = 0° 
is  independent  of  the  y-coordinate : 


2 -ikx  sin  0 
™ Tra  o 

p2-3?Ie 


-ikz  cos0. 


COS0, 


i27r(x/d^)  -kz^ 

[(  27t /kd  L ) 

»-sineo] 

I2'1] 

1/2 

* 

[(  27r/kd  x) 

- sin0  1 
oJ 

2 - 

\ 


i 


(27) 


Consider  the  case  kz  = 0,  then 


10 


WP22-1-42011 


lP2l  " 

1 

i 

e 

27r(x/d1)+7r/2 

COS0„  f 

0 1 

(27r/kd^)  - sin0Q 

2 

(28) 


The  two  terms  in  Eq.  (28)  can  be  considered  as  vec- 
tors in  the  complex  plane . Thus  we  will  have  a maximum  when  the 
phase  2v  (x/d^)  + v/ 2 of  the  second  term  becomes  zero,  or  when 
the  x-coordinate  is 


V 

max 


-dL/4, 


(29) 


independent  of  9o,  d2>  or  a. 

Similarly  a minimum  will  be  obtained  at 

xmin  = + V4- 


(30) 


The  positions  of  the  minima  and  the  maxima  agree  well 
with  the  data  in  Fig.  3 for  0q  >0°.  For  the  values  at  the  maxi- 
mum and  at  the  minimum  Eq . (28)  yields  1.19  and  0.21  at  0 = 60°; 
and  2.74  and  1.29  at  9 = 80°,  which  are  in  a reasonable  agree- 
ment with  the  values  shown  in  Fig.  3. 

Eq.  (28)  yields  the  maximum  pressure 
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under  the  condition  kz  = 0.  However,  for  kz  > 0 the  second  term 
in  the  more  general  Eq . (27)  is  exponentially  decreased;  there- 
fore, even  if  we  consider  all  values  of  x,y,  and  z,  the  absolute 
maximum  of  | ?2|  must  be  on  the  plane  kz  = 0,  and  its  value  is 
given  by  Eq . (31),  when  0Q  = 0°  , > 0°  , kd^  < tt,  and  kd2  < kd^. 

Fig.  4 shows  |p|  vs.  kx  for  ka  = 0.5,  1.0,  and  1.5, 
and  kd^  = kd2  =3,  ky  = kz  = 0,  0Q  = 0° , 9q  = 60°.  The  pressure 
levels  increase  as  ka  increases  because  a larger  fraction  of  the 
array  plane  radiates.  A larger  fraction  of  the  array  plane  radi- 
ates also  if  ka  is  kept  constant  and  kd-^  and  kd2  are  decreased. 

| P J vs.  kx  for  the  latter  case  is  shown  in  Fig.  5,  where  ka  = 1, 
kdL  = kd2  = 2,  2.5,  and  3,  and  kz  = ky  = 0,  0Q  = 0°  0Q  = 60°. 

Note  that  for  kd^  = kd2  < 3 in  Fig.  5 more  than  one  cycle  in  the 
undulatory  pressure  is  shown,  because  due  to  symmetry  the  pressure 
values  are  repeated  at  x,  x + d^,  x + 2d^,  etc.  The  pressure  dis- 
tribution becomes  relatively  more  undulatory  as  the  spacing  bet- 
ween pistons  is  increased.  When  ka  = 1.5  and  kd^  = kd2  = 3 in 
Fig.  4,  and  when  ka  = 1 and  kd^  = kd2  = 2 in  Fig.  5,  then  the 
pistons  a^e  touching  each  other,  and  the  same  fraction  tt/4  of 
the  array  plane  radiates  in  both  cases.  However,  the  maxi- 
mum pressure  in  Fig.  4 for  the  larger  piston  ka  = 1.5  is 
about  24%  larger  than  for  the  smaller  piston  ka  ^ 1 in  Fig.  5. 
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The  pressure  is  also  more  undulatory  (i.e.,  the  ratio  of  the 
maximum  pressure  to  the  minimum  pressure  increases)  for  the 
larger  piston.  All  of  the  above  features  can  be  explained  quali- 
tatively by  considering  the  two-term  approximation  P ^ As  kd^ 
is  increased,  the  magnitude  of  the  second  term  in  the  parentheses 
in  Eq  • (28)  is  increased;  therefore  the  pressure  be- 

comes more  undulatory  as  the  spacing  between  pistons  is  increased. 

2 

If  we  keep  the  ratio  ttbl  /d^2  constant  in  Eq . (28)  or  Eq . (31), 
then  for  a larger  piston  we  have  a larger  kd-^,  and  thus  a larger 
maximum  pressure  and  a more  undulatory  pressure  distribution. 

Although  ?2  is  useful  as  a guide  to  the  general  be- 
havior of  the  pressure,  it  is  of  less  value  in  the  calculation  of 
the  actual  pressure  magnitudes;  see  Table  I.  Furthermore,  ?2 
does  not  agree  with  the  behavior  of  the  pressure  for  ka  - 0.5 
in  Fig.  4;  while  the  maxima  and  minima  should  be  located  at 
+d^/4  independent  of  0o>  d2,  or  a,  the  maxima  and  minima  have 
been  shifted  for  the  ka  = 0.5  curve  in  Fig.  4.  This  curve  is 
the  only  one  of  those  shown  in  Figs.  3,  4,  and  5,  for  which 
a < d^/4,  i.e.,  for  which  Eq . (29)  predicts  a maximum  pressure 
on  the  adjacent  baffle  off  the  surface  of  the  piston;  however, 
the  actual  maximum  is  on  the  surface  of  the  piston. 

Moreover,  according  to  Eq . (27)  P2  is  independent  of 
ky,  while  the  contours  of  |P|  in  Fig.  6 for  z = 0 show  a ky  de- 
pendence. The  pressure  contours  in  the  planes  x = 0 and  y = 0 
are  shown  in  Figs.  7 and  8,  when  kd^  = kd£  = 3 , ka  = 1 , 0n  = 0° , 
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and  0Q  = 60° . Since  0Q  = 0° , the  pressure  field  is  symmetric 
about  the  x-axis,  and  only  pressures  for  ky  > 0 are  shown.  For 
large  values  of  kz,  |p|  = 0.63  in  Figs.  7 and  8. 

C . RADIATION  IMPEDANCE 

th 

The  force  on  the  m ,n  piston  divided  by  the  velo- 

o 

city  v , is  the  radiation  impedance  pcra  Z,  where  Z is  a dim- 
ensionless radiation  impedance  coefficient.  Due  to  symmetry  Z 
is  the  same  for  all  pistons. 

9 2tt  a 

Z = (l/7ra  ) r da  | dr  r P = R + iX 

Jo  Jo  z=0  (32) 


where  x = r cosa,  y = r sina,  and  P is  given  by  Eq . (21).  The 
integrals  in  Eq . (32)  are  the  same  as  in  Eqs.  (19)  and  (20),  thus 


^ w(n,m)ka 
w(n,m)ka 


2 


1 

T(n,m) 


V 


J 


(33) 


For  kd^  < it,  kd2  < tt  the  dimensionless  radiation  resistance  R 
is  given  by  the  n = m = 0 term. 


2j^(ka  sin0Q)  ^ ^ 

ka  sin0„  cos6 


(34) 
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while  the  reactance  X is  given  by  the  infinite  sum  without  the 
n = m = 0 term.  The  corresponding  quantities  zc  = Rc  + 1XC  for 
the  continuous  plane  radiator  are^ 

Rc  = l/cos0o,  (35) 

Xc  = 0.  (36) 


I 

I 

I 

I 

I 

I 

I 


t 


f 


The  fact  that  the  radiation  resistance  should  be  proportional 
to  l/cos©o  can  also  be  deduced  for  a finite  array  for  0Q  < tt/2 
by  the  use  of  arguments  involving  the  array  beamwidth  in  the 
far  field.® 

One  could  also  obtain  the  radiation  impedance  by 

summing  an  infinite  series  of  the  appropriate  mutual  coupling 
19 

coefficients  ; however,  the  series  in  Eq . (33)  converges  more 
rapidly  than  such  a series  of  mutual  coupling  coefficients,  and 
thus  is  more  useful.  An  approximate  expression  for  Z has  been 
obtained  beforezu  by  the  application  of  Poisson'  s sum  formula 
to  a summation  of  approximate  mutual  coupling  coefficients,  but 
the  exact  Eq.  (33)  is  even  simpler  than  the  approximate  expression. 

Fig.  9 shows  X vs.  kd^  = kd2  for  three  values  of  ©0, 
and  ka  = 0.5,  0O  = 0°.  As  the  spacing  between  pistons  is  in- 
creased, the  radiation  resistance  R is  decreased,  which  is  ob- 
vious from  Eq.  (34),  and  the  radiation  reactance  X is  increased, 
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as  shown  in  Fig.  9.  The  radiation  resistance  is  proportional 

to  the  energy  which  is  radiated  to  the  far  field;  the  radiation 

reactance  is  proportional  to  the  energy  which  just  fluctuates 

21 

in  the  near  field.  When  the  whole  plane  is  vibrating,  then 
there  is  no  energy  fluctuating  in  the  near  field,  see  Eq.  (36), 
while  for  the  set  of  parameters  used  in  Fig.  9 the  near  field 
fluctuating  energy  is  increased  as  the  ratio  of  the  active  pis- 
ton areas  to  the  passive  baffle  area  is  decreased  by  increasing 
kd-^  = kd2  and  keeping  ka  constant. 

If  the  ratio  of  the  active  piston  areas  to  the  passive 
baffle  area  is  decreased  by  decreasing  ka  and  keeping  kd^  = kd2 
constant,  the  near  field  fluctuating  energy  may  be  increased  or 
decreased,  see  Fig.  10,  where  X vs . ka  is  shown  for  three  values 
of  9q,  and  kd-^  = kd2  = 3,  0Q  = 0°  . 

Fig.  11  shows  X vs . 0Q  for  three  values  of  ka,  and 
kdx  = kd2  = 3,  0Q  = 0°. 

In  Fig.  12  the  radiation  impedance  of  a piston  in  the 

infinite  array  is  compared  with  the  average  radiation  impedance 

R + iX„  of  a piston  in  a finite  array  of  12  rows  (parallel 

to  the  y-axis)  and  229  columns  (parallel  to  the  x-axis) . The 

radiation  impedance  of  a piston  in  the  finite  array  is  obtained 

19 

by  summing  the  appropriate  mutual  coupling  coefficients,  and 
in  general  the  impedances  vary  from  piston  to  piston.  The 
average  radiation  impedance  for  the  finite  array  is  obtained 
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by  summing  the  impedances  of  all  the  pistons  and  by  dividing  the 
sum  by  the  total  number  of  pistons.  The  agreement  in  Fig.  12 
for  kd^  = kd2  = 2.51,  ka  = 0.80,  0Q  = 0°,  is  good  except  near 
0Q  = 90°  where,  of  course,  R for  the  infinite  array  becomes  in- 
finite, while  for  the  finite  array  Rflv  is  finite. 


D.  CONCLUSIONS 


We  have  obtained  more  rapidly  convergent  infinite  series 
expressions  for  the  near  field  pressures  and  the  radiation  im- 
pedances of  circular  pistons  in  an  infinite  phased  array.  Numer- 
ical results  have  been  obtained  for  a limited  number  of  cases, 
and  these  cases  indicate: 

1)  A single-term  approximation  P^,  see  Eq . (24),  is 
sufficient  to  represent  the  pressure  if  the  distance  from  the 
array  kz  is  larger  than  1 or  if  the  array  consists  of  small 
(ka  « 1)  pistons  close  together. 

2)  A two-term  approximation  P2,  see  Eq . (27),  gives 
a rough  estimate  of  the  pressure  when  0Q  = 0°,  0Q  >0°,  kd^  < tt, 
and  kd2  < kd^;  in  particular  an  approximate  location,  see  Eq . (29), 
and  value,  see  Eq.  (31),  of  the  maximum  near  field  pressure  can 

be  obtained,  and  some  gross  features  of  the  pressure  distribution 
can  be  predicted. 

3)  For  0Q  = 0°  the  pressure  distribution  is  relatively 
smooth,  but  as  0Q  increases,  the  pressure  distribution  becomes 
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more  undulatory,  and  the  maximum  pressure  increases. 

4)  If  we  keep  the  piston  radius  constant  and  in- 
crease the  spacing  between  pistons,  the  pressure  distribution 
becomes  more  undulatory. 

5)  If  we  keep  the  ratio  of  the  piston  areas  to  the 
baffle  area  constant  (sometimes  called  the  "packing  factor") , 
then  larger  pistons  further  apart  have  a larger  maximum  pressure 
and  a more  undulatory  pressure  distribution  than  smaller  pistons 
closer  together. 

6)  The  near  field  pressure  magnitude  and  the  radi- 
ation resistance  increase  without  bounds  as  9q  ->  90°  , while  the 
radiation  reactance  remains  finite  for  kd^  < tt  and  kd2  < v. 

7)  The  radiation  impedance  of  a piston  in  the  in- 
finite array  agrees  well  with  the  average  radiation  impedance 
of  a piston  in  a large  finite  array,  except  near  9q  = 90°. 
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TABLE  I 

Comparison  of  | P2 j max  and  the  actual 
maximum  of  ) p|  . 


kdL 

kd2 

ka 

eo 

urn 

Maximum  | P | 

3 

3 

mm 

60° 

1.19 

1.18 

3 

3 

■■ 

00 

0 

e 

2.74 

2.55 

2.5 

2.5 

■■ 

60° 

1.39 

1.26 

2 

2 

■■ 

60° 

1.95 

1.64 

3 

3 

H 

60° 

2.67 

2.04 

3 

3 

0.5 

60° 

0.30 

0.47 

3 

2 

1 

0 

0 

v£> 

1.78 

1.59 

0=0°  for  all  cases 
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SEVERAL  ka. 


Ky  * kz  * 0 


FIGURE  4.  P vs.  kx  FOR  SEVERAL  ka 
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FIGURE  5.  P vs.  kx  FOR  SEVERAL  kd 
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FIGURE  8.  PRESSURE  CONTOURS  IN  THE  PLANE  y 


FIGURE  9.  X vs.  kd,  = kd-  FOR  SEVERAL  0 
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FIGURE  10.  X vs.  ka  FOR  SEVERAL  6 
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FIGURE  II.  X vs,  FOR  SEVERAL  ka 


